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ABSTRACT 

Inspired by the recent identification of the first candidate BH-WD X-ray binaries, where the compact accre- 
tors may be stellar-mass black hole candidates in extragalactic globular clusters, we explore how such binaries 
could be formed in a dynamical environment. We provide analyses of the formation rates via well known 
formation channels like binary exchange and physical collisions and propose that the only possibility to form 
BH-WD binaries is via coupling these usual formation channels with subsequent hardening and/or triple for- 
mation. Indeed, we find that the most important mechanism to make a BH-WD X-ray binary from an initially 
dynamically formed BH-WD binary is triple induced mass transfer via the Kozai mechanism. Even using the 
most optimistic estimates for the formation rates, we cannot match the observationally inferred production 
rates if black holes undergo significant evaporation from the cluster or form a completely detached subcluster 
of black holes. We estimate that at least 1% of all formed black holes, or presumably 10% of the black holes 
present in the core now, must be involved in interactions with the rest of the core stellar population. 
Subject headings: X-rays: binaries — star clusters: general — stars: kinematics and dynamics. 



1. INTRODUCTION 

The question of whether black holes (BHs) are present in 
globular cluste rs (GCs) has been discussed extensively in the 
literat ure (e.g..lKulkarni et alJl!993b ISigurdsson & Hernauisl 
19931: iPortegies Zwart & McMillanlboOOt iMiller & Hamilton 
2002bl) . The most obvious way to detect a BH in a GC 
is if the BH is in an X-ray binary. IKalogera et al.l (12004b 
showed that, if a BH X-ray binary is formed dynamically 
by an exchange interaction in the core of a dense globu- 
lar cluster, it is very unlikely to be detected, as the duty 
cycle for such binaries is extremely low. Tidally captured 
BH binaries, in contrast, would be continuously luminous, 
and the lack of Galactic BH X-ray binaries implies that any 
such tidal captures destroy the potential companions. We 
note that only non-d egenerate donors were considered by 
IKalogera et al.l (12004 . Indeed, in Galactic GCs, where the 
total number of LMXBs is relative ly small, no BH X-ray b i- 
nary has been found so far (e.g.. IVerbunt & Lewinl 120061) '). 
Among low-mass X-ray binaries (LMXBs) in the globular 
clusters of early-type galaxies, however, likely BH X-ray bi- 
nary candidates h ave b een reported (e.g.. lAngelini et al.l2doTt 
Pi Stefano et al Il2002t [Kundu et al.l2002HSarazin et al.l2003l: 
Kim et al 1 120061) . The strongest such candidates have X-ray 
luminosities Lx > 10 39 ergss" 1 , well above the Eddington 
limit for a neutron star accreting helium, and are known as 
ultraluminous X-ray sources (ULXs). 

A particularly interesting ULX has been identified in a 
globular cluster in NGC 4472. It has L x ~4x 10 39 ergs s" 1 
and has shown strong variability (likely due to variable ab- 
sorption), indicating that it is a single source and thus likely a 
black hole X-ray binary dMaccarone et al. 2007). Keck spec- 
troscopy of the globular cluster (RZ 2109) associated with 
this source identified strong, broad (2000 km/s) [O III] emis- 
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sion lines, interpreted as an outflow from a stellar mass black 
hole accreting above or near its Eddington limit dZepf et al.l 
120081) . The low Ha/[0 III] ratio suggests a hydrogen-poor 
(white dwarf) donor, while the assumption that the breadth of 
the line is due to a wind driven fr om near-Eddington accretion 
implies a BH mass of 5 - 20M Q dGnedin et al.ll2009l) . 

The ultraluminous X-ray source CXO J03383 1.8-352604 
is associated with a globular cluster in the Fornax elliptical 
galax y NGC 1399, with L x > 2 x 10 39 ergs s _1 dlrwin et al.l 
2009). This source shows strong [O III] emission lines (less 
broad; a ~ 70 km/s) and little or no hydrogen emission. It 
may be another black hole- white dwarf X-ray binary accreting 
near Eddington, although llrwin et all d2009) suggest a tidal 
disruption of a white dwarf by an intermediate mass black 
hole (IMBH). 

How common are such systems? iKim et al.l (120061 K06) 
found 8 such ULXs, with L x > 10 39 ergs s _1 , in 6173 glob- 
ular clusters, where these GCs have an average mass of 
6 x 1O 5 M . Humphrey & Buotd (120081 HB08) found 2 such 
objects among 3782 globulars with total Ly ~ 10 9 L Q (to- 
tal mass M~ 3x 10 9 Mq). These surveys suggest that GC 
ULXs are present at the rate of 2.0+1% x 10" 9 and 7^ 5 x 10" 10 
per M q, respectively , using Gehrels' statistics at 90% confi- 
dence (Gehrels 1986). From the overlap in these numbers, we 
estimate that there is ~1 such object per 1 - 2 x 10 9 M Q in 
GCs. 

These surveys were not completely independent, as sev- 
eral galaxies (including NGC 1399) were contained in both 
surveys. Some fraction of the observed globular ULXs may 
not be BH-WD systems, but accreting IMBHs. The simi- 
larities between t he two well-observed ULXs, discussed in 
llrwin et al.l (120091) . suggest that they are either both BH-WD 
binaries or both accreting IMBHs. For this paper, we explore 
the consequences of the assumption that all globular cluster 
ULXs are BH-WD binaries, in particular whether such num- 
bers can be produced, and by what mechanisms. 

Let us consider now what this frequency implies for the 
formation rate of such systems per BH. In globular clusters, 
roughly each 150-200 M Q of currently remaining stellar mass 
produced a BH in the past, where half of these BHs have 
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masses above 10M Q (here, less cluster mass per BH is re- 
quired in a metal-poor cluster). The retention fraction of 
BHs immediately after their formation, if supernova kicks are 
taken into acco unt, is about 30-40%, for the escape veloc- 
ity of 50 km/s (Be lczvnski et al.ll2006l) . We estimate there- 
fore that, if no dynamical interacti ons between stars occur, a 
present day average cluster in the iKim et all (120061) sample 
(6 x 10 5 M Q ) would retain about 500-800 BHs with masses 
above 1OM . 

Standard consideration of the consequences of dynamical 
interactions between BHs and other stars in clusters assumes 
that there is a tendency towards equipartition, where BHs 
quickly detach from other stars and for m a BH subclu ster (this 
is also known as the Spitzer instability, ISpitzerll969l) . Within 
this subcluster they interact only with each other and, as a re- 
sult, fairly quickly e vaporate from the clu ster. Assuming the 
evapor ation times of Wa tters et al.l (120001) and lKalogera et al] 
d2004l) . only one BH could remain per cluster at the current 
epoch. This evaporation timescale assumed simple evapora- 
tion from an isolated cluster and is probably on the high side. 
Detailed numerical calculations of BH subcluster dynamical 
evolution, however, have shown that in some circumstance a 
significant fraction of BHs will remain at 10 Gyr time - in 
massive clusters, the fraction of BHs remaining in the cluster 
at 10 Gyr is ~ 20% o f the BHs retained after their formation 
dO'Learv et al.l [2006). It has been also noticed that in mas- 
sive clusters the BH subcluster does not reach equipartition 
at Hubble time. Sim ilar results were found in other studies 
dDowning et al.l2009h . where a whole globular cluster, includ- 
ing normal stars, was modeled using a Monte Carlo treatment. 
Even though black holes were found to be strongly centrally 
concentrated and interacting mainly with each other, a signif- 
icant fraction of initial BHs (up to 25%) not only remained in 
the cluster but also participated in interaction with other stars. 
(It is not clear which fraction is purely in the BH subclus- 
ter and which fraction did not detach from the rest of stellar 
population.) We assume therefore that the number of massive 
BHs (> 10M Q ) retained, and available for intera ctions with 
core st ars, could be from one per average cluster in lKim et alj 
(120061) sample (i.e. roughly 4000-6000 BHs in each whole 
observed sample, K06 and HB08) up to 125-200 per average 
K06 cluster or 160-260 per average HB08 cluster (i.e. about 
5 x 10 s - 1.2 x 10 6 BHs in the observed samples). The up- 
per limit corresponds to the case when the number of massive 
BHs is 10% of all initally formed massive BHs. 

The lifetime of the persistent mass transfer (MT) at the Ed- 
dington level (~4x 10 39 ergs s" 1 ) from a WD to such a BH, 
once started, is tx « 2 x 10 5 yr and is about 3 times longer 
for ULX luminosities Lx > 10 39 ergs s . If all but one BH 
is evaporated, the formation rate of BH-WD binaries must be 
extremely high - up to one X-ray binary formation per BH 
per Gyr. The minimum required formation rate, when a good 
fraction of BHs is retained which also do not form a BH sub- 
cluster, is about 4 x 10~ 3 per BH per Gyr (and more likely 
10~ 2 per BH per Gyr, if many clusters are metal -rich and have 
lower escape velocity than 50 km/s). We therefore explore 
which dynamical formation channel would be able to provide 
the formation of BH-WD X-ray binaries at rates of 4 x 10~ 3 
to 1 per BH per Gyr in an average dense massive cluster. We 
can also define this as 4 x 10~ 3 //bh,o.i per BH per Gyr, where 
/bh.o.i = 0.1/BH,tot is the fraction of BHs that is retained in the 
cluster and not detached in the BH subcluster, and is normal- 
ized to 10% of all initially formed massive BHs (> 10M Q ). A 



smaller rate is possible only if essentially all globular clusters 
have experienced significant tidal stripping and the number of 
remaining BHs per present stellar mass unit is significantly 
higher than 1 BH per 150 M@. To analyze the dynamical 
formation, we will proceed in reverse order: first, we will 
consider which BH-WD binaries, once formed, can become 
X-ray binaries, and then we will consider at what rates these 
BH-WD binaries can be formed via different dynamical chan- 
nels. 

2. FATE OF A BH-WD BINARY 

2.1. Characteristic times 

In a dense stellar system, once a BH-WD binary is formed 
through a dynamical encounter, its further binary evolution 
could be affected by subsequent encounters with other stars. 
The cross section for an encounter between two objects of 
total mass m tol with a distance of closest approach less than 
'"max is computed as 
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where the second term accounts for gravitational focusing, 
with v 2 = 2Gm tot /r max and Voo the relative velocity at infin- 
ity. For strong interactions, r max is usually on the order of the 
binary semimajor axi s: r max = ka, where k is of order unity 
dHut & Bahcalllll983l) . In the limit of strong gravitational fo- 
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In our case, the first object is the BH-WD binary of mass 
twbhwd = wbh+mwd while the second object is a core star of 
mass m*. In globular clusters, the close approaches that we are 
interested in have v 2 ^> v^, and the BH mass is significantly 
more massive than both its WD companion and a typical core 
star. Thus, the rate at which a BH-WD binary undergoes a 
strong (binary-single) encounter is 
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The final equal sign in equation ([3]i defines the dimensionless 
parameter K. The time-scale for a BH-WD binary to expe- 
rience a strong encounter can be calculated as tbs = 1 /Tbs = 
10. 6 K _1 Rq /a Gyr: see FigureQ] Throughout this paper, we 
consider clusters with core number densities n c near 10 5 pc~ 3 , 
velocity dispersions of ~ 10 km s" 1 , and black hole masses of 
~ 15M Q : consequently, K is of order unity. 

Once formed, the orbit of a BH-WD binary starts to shrink 
due to gravitational radiation. The time r gw for the orbit to 
decay enough for mass transfer to commence depends on the 
initial po st-encounter b inary separation a and eccentricity e. 
Applying Pejeri ( 119641) equations, we can find the maximum 
semimajor axis a, as a function of post-encounter eccentricity 
e, for which a system will start MT within any specified time, 
e.g., 1 Gyr and 10 Gyr (see the dotted curves in FigQ]). 
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companion exchange, or a physical collision. Such outcomes, 
which we now consider, could potentially prevent a BH-WD 
binary from ever reaching the MT stage. We note that if an en- 
counter occurs with another binary, a possible outcome is also 
a hierarchically stable triple system; we consider this option 
in the next subsection. 

The collision cross section for binary-single interactions is 
approximately (for a non-eccentric binary) 5 : 
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where V} = Gm tot /(2a) and the maximum radius R among the 
encounter participants is pr esumed to be less th an a. Our def- 
inition of V g follows that of lHeggie et alj d 19961) . We note that 
V g = v c in the case of three equal mass stars. We have veri- 
fied with numerical experiments for many mass combinations 
that the use of V} yields the appropriate scaling of the cross 
section with mass (provided Voo <C Vg). 

Using equations © and (0, we can now estimate the frac- 
tion of binary-single encounters that result in a collision: 



FIG. 1. — Orbital parameters (initial separation ci/Rq and eccentricity e) 
of a 15 Mq BH and 0.6 Mq WD binary system. Lines of constant merger 
time T gw (dotted curves) and encounter time rgs with a single star (dashed 
lines) are shown at specified values. Also shown is the separation curve o SC p 
for which r gw = res (solid curve), assuming the dimensionless parameter K 
defined by equation (5) equals 1 . 

In Figure Q] we also show the semimajor axis a sep (e) such 
that tbs = T gw (e), assuming K = 1 . BH-WD binaries with post- 
encounter semimajor axes a{e) < a se p(<?) will become BH-WD 
X-ray binaries. BH-WD binaries with post-encounter semi- 
major axes a(e) > a sep (e) will have one or more strong en- 
counters before they shrink to the point of merger due to grav- 
itational waves radiation. In the next subsection, we consider 
the possible outcomes of such encounters. 

2.2. Encounters with single stars 

Here we consider the outcomes and their frequencies for en- 
counters between a BH-WD and a single star. An encounter 
between a soft binary and a third star can lead to ionization 
if the third star approaches with sufficient speed; specifically, 
for ionization to be energetically possible, the relative veloc- 
ity at infinity between the binary and the third star must 
exceed the binary's critical velocity v c defined such that the 
total energy of the binary-single system is zero: 

2 m BHWD + »J* %D CfflfiH 



Ccoll 

a 



1 (215R 
4k [ a 



(6) 



'WBHWD 



(4) 



For mBH = 15Mq, mwD = 0.6M©, and m* = Q.6Mq, the critical 
velocity v c is less than 10 km s -1 only if the semimajor axis 
a of the binary exceeds ~ 3 x lO 4 /?©. Even for encounters 
between a BH-WD binary and a 15M© black hole, the critical 
velocity v c is less than 10 km s" 1 only for a > 2000,/?©. Con- 
sequently, in a typical cluster with velocity dispersion ~ 10km 
s -1 , ionization by itself can not destroy BH-WD binaries that 
are able to start MT within 10 Gyr, as all such binaries are 
easily of sufficiently small semimajor axis (see Fig. [TJ. Es- 
sentially, a BH-WD binary with a < 2000/?© can be conser- 
vatively classified as a hard binary. 

A strong encounter between a hard binary and a third star 
would ultimately lead to preservation of the BH-WD binary, a 



where k is the ratio of the maximum pericenter considered as 
an encounter to the semimajor axis a of the binary. Equation 
(O implies that if a non-eccentric BH-WD has a < 5Rq, then 
all encounters with R = 0.6/?© MS star at k < 2 will result in 
a merger, while for encounters with k < 1, a < 15/?© will re- 
sult in a merger. Consequently, we expect that most BH-WD 
binaries that are able to start MT within 10 Gyr (see Fig. Q~|i 
will experience a merger if a binary-single encounter with a 
MS star occurs. 

To examine the collision cross s ection more carefully , we 
perform scattering experiments as in lFregeau et al.l ( 120041) . but 
for our BH-WD binary with the separations from 1 to 100 
Rq and fixed v,^ = 10 km/s. We find the fc co ii, defined by 
equation (f2]i with a = cr co n and k = k m \\, varies from ~ 1 for 
a = IRq to ^0.1 for a = 100/?©. For example, £ co u w 0.24 
for a = IORq, and, accordingly from equation (0, the merger 
rate for such BH-WD binaries that are able to start MT within 
10 Gyr exceeds 0.2 per Gyr. We note that effectively this rate 
could be a factor of few smaller as a fraction of collisions will 
result in a merger between a BH and a MS star. We have 
verifed with hydrodynamicsl simulations that at least some 
of these collisions will not affect strongly the initial BH-WD 
binary. 

The relative fraction of WDs in the core, at the age of 
several to 14 Gyr, is /wd ~ 0.2 of the total core popula- 
tion, and Tbs.wd = /wdTbs- We also note that some simu- 
lations have shown that WDs can contribute up to 70% of 
the total co re population, as was found in simulations per- 
formed for iFregeau et all d2009bl) without WD birth kicks. 
From eq. ©, we roughly estimate that a physical collision 
with a R = 0.01/?© WD would occur only in binaries with 
a < O.IRq, which is well below a sep (e). A more likely out- 
come in such cases would be preservation or a companion ex- 
change. Exchanges preferentially occur if the incoming star is 
more massive than pre-encounter companion and accordingly 
could happen if BH-WD binary had a light WD companion. 

5 This approximation was derived using results from Fregeau et al. 1 2004). 
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FIG. 2. — Fraction of encounters, for k = 20, that result in a triple forma- 
tion. Encounters are between a BH-WD 15Mq +0.6Mq binary with sep- 
arations a i = 15,35, and 80Rq (green triangles, red stars and blue squares 
accordingly) and 0.6Mq + 0.6Mq binary with separations «2- Each point 
shows the result for 10 000 encounters. The solid line is the fitting with 

0.29(l+e ( - 31 * !(n 2/' , i )+1, ) _1 

As a result of exchange, the post-encounter binary separation 
will be increased by the ratio of new companion mass to the 
old companion mass. To reiterate, we consider here only the 
cases when a BH-WD binary has a(e) > a sep (e). Therefore, 
the resulting binary will most likely never be able to reach 
MT, as its binary separation will exceed the value necessary 
for MT to start within 10 Gyr. If a binary is preserved, though, 
it will most likely then experience a consequent encounter 
with a single star. Even if WDs make up as much as 70% of 
the core population, this consequent encounter has high prob- 
ability to have a MS star as a participant and therefore to result 
in a merger. We conclude that strong encounters with single 
stars will not create BH-WD binaries with a(e) < a S ep(<?) and 
therefore will not lead to a direct formation of a BH-WD X- 
ray binary. 

2.3. Encounter with binary stars: role of triples 

If a BH-WD has an encounter with another binary, a possi- 
ble outcome is a hierarchically stable triple. For two identical 
binaries, the maximum impact parameter for which strong en- 
counters st ill could occur i s a bit larger than for a binary-single 
encounter (Mikkola 1983): 



ctbb ~4.67r(ai+a 2 ) 
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Here v c is the critical velocity of the two binaries with the 
masses m\ and ni2'. 
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where the reduced mass ij, = m\m2/(m\+m2), where /nn, mn, 
m.2i and m.22 are the masses of the binary components, and 
where a\ and «2 are the semimajor axes. 

For hard binaries with equal masses and binary separations, 
~ 25% of all binary -binary encounters within maximum ap- 



proa ch described by eq. Q and vl/v 2 ^ = 10 form a triple sys- 
tem (lMikkolalll983t IFregeau et al.l 120041) . with triple forma- 
tion increasing with v^/v^. For a larger mass ratio between 
binaries, or for a softer incoming binary of less mass, this 
fraction can be even higher. 

We performed numerical experiments in which BH-WD bi- 
naries consisting of a \5Mq BH and a 0.6M© WD encoun- 
tered at Vqo = lOkm/s binaries consisting of 0.6 M© compan- 
ions. We record the fraction of encounters that result in a 
stable triple with an inner binary consisting of the initial BH 
and a companion 0.6M© star. We run the simulations for 3 
values of BH-WD binary separation and 11-17 values of the 
incoming binary separations, performing 10000 encounters 
for each. For the nu merical integration, we use Fewbody 
( IFregeau et alj [20041) . a small A^-body integrator. The re- 
sults demonstrate that periastrons somewhat larger than those 
suggested by eq.© can result in a triple formation, and we 
choose the maximum possible periastron for the encounters 
to be fairly large, r max p = 20(ai +02) (that is, k = 20). The im- 
pact parameters for every encounter are distributed uniformly 
in area. The eccentricities are distributed thermally. The re- 
sults are shown at the Figj2] 

We find, that for almost all the encounters in which the in- 
coming binary is wider than the BH-WD binary, the efficiency 
of triple formation is ~ 20-30% with k = 20. This gives 
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The formation rate of triple systems, for hard binaries, com- 
pared to to binary-single encounters rate ([3]) with k = 2, is then 
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Here / w b is the fraction of stars that are hard binaries, but 
still wider than the BH-WD binary. In deriving the above 
equation, we have applied equipartition of energy, so that the 
Voo of a typical binary in the core is y/2 times less than that of 
a typical single star. 

The overall binary fraction in observed dense, but not core- 
collapsed, clusters at the current epoch, is about a few per cen t 
dAlbrow et al.l200Ulrvanova et alj2005atlMilone et al.l2008l) . 
With / w b =5% and a flat distribution of binary separations 
between, e.g., 20 and 20007?©, a BH-WD binary that has a = 
2QRq can form a triple about 30 times per Gyr. This exceeds 
the interaction rate with single stars. 

Although these triples are stable in isolation, they can be 
destroyed during their next dynamical encounter. Before this 
next encounter occurs, the presence of the third component 
can significantly affect the eccent ricity of the inner BH-WD 
binary via the Kozai mechanism dKozaill 19621) . provided the 
inclination of the outer orbit is large enough (sin;'o > (2/5) 1 / 2 , 
or equivalently ;'o > 39°, the Kozai angle). The maximum 
eccentri city that can be achieved via the Kozai mechanism 
is (e. g., Ilnnanen et alJ 11997b lEggleton & Kiseleva-Eggletonl 
12001 
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In dynamically formed triples, the inclinations are dis- 
tributed almost uniformly, and so the number of Koz ai triples 
is expected to be proportional to the solid angle (e.g.. lIvanoval 
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l2008ai where, depending on the sample of runs, the triples af- 
fected by the Kozai mechanism constituted from 30% to 40% 
of all formed triples, while a uniform distribution predicts 
37%). Thus, the fraction of the dynamically formed triples 
with a post-encounter inclination greater than ; is f t = 1 — sine. 

The BH-WD binaries that could become MT systems are 
only those that have a less than the maximum a sep (< 80/?© 
for e < 0.99). To find the formation rate of such triples, 
we first determine the minimum value of e max such that a 
triple with the specific binary separation a must achieve: 
this is found by inverting the numerically obtained a sep (e) at 
k = 2. The obtained required e ma x then provides the fraction 
of formed triples that can be bring the inner binary into MT 
via the Kozai mechanism: 



f . = l-sin^l-,/3/5^+2/5. 



(12) 



We define triple induced mass transfer (TIMT) systems to 
be those binaries that are brought into MT via the the Kozai 
meachanism. We find that the fraction of TIMT systems 
ranges between 0.07 (fora= 15/? Q ) and 0.01 (for most separa- 
tions, a > 4QRq). Combining it with the triple formation rate, 
we find that all BH-WD binaries with 15R Q < a x < 80R Q 
could make it to the MT through the triple formation mech- 
anism. At a conservative level, assuming even that all the 
binaries that have a binding energy 10 times more than the 
kinetic energy of an average object in the core are destroyed, 
we find that a BH-WD binary with a x > 50R Q has a 10-40% 
chance to become become a TIMT system within 1 Gyr (10% 
for «i ~ 80R Q and 40% for a x ~ 50R & ); this chance is 100% 
for «i < 35Rq. During several Gyr, all BH-WD binaries with 
«i < 80^0 can become a TIMT system at least once. 

This TIMT system formation will be successful only if the 
time between subsequent encounters is longer than the time 
necessa ry to achieve e max . The period of the cycle to a chieve 
e max is dlnnanen et ai]|1997tlMiller & Hamirtonll2002ah 
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where m ;l and ma are the companion mases of the inner bi- 
nary, e\ is the initial eccentricity of the inner binary, m is the 
mass of the outer star, a is the initial semimajor axis for the 
outer orbit, and b D = a a (l -e 2 ) 1/2 is the semiminor axis of the 
outer orbit. Roughly, a will be of the order of magnitude 
of «2- Considering the extreme case a 2 3> a.\, for which r Koz 
would be maximum, we find: 
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Here f b is the binary fraction, r BB = T BB = (o- B B./b«cVoo) l , 
and r B s = 1 /T B s has been evaluated at k = 2. 



An av e rage dy namically formed triple would have e a w 0.9 
llvanoval d2008bl) . For the purpose of the estimate, we adopt 
that a BH-WD binary would have its initial eccentricity dis- 
tributed thermally, with an average e\ ~ 2/3 and adopt ;'o = 54° 
(this is the mean inclination of dynamically formed Kozai af- 
fected triples). From eq.dToll. we have then that triples with 
15#© < fli < 80/?o and a 2 < 10 5 R e (where 4500/?© is the 
boundary between hard and soft binaries) would have their 
Kozai time significantly shorter than their collision time with 
either binary or single stars. We conclude that once a potential 
TIMT system is formed, it will succeed in bringing its inner 
BH-WD to the mass transfer before its next encounter. Ac- 
cordingly, all BH-WD binaries with a < 80Rq will become 
X-ray binaries via TIMT mechanism. 

2.4. Multiple encounters: hardening 

It is well known that thro ugh multip le fly-by encounters 
hard binaries get harder (e.g.. lHuj|1983l) . The number of en- 
counters necessary to change th e binary energy, e.g., by a fac- 
tor of four, can be estimated as (Heggie & Hut 2003) 



AWi«log(4)/log(l+<5), 



(16) 



where S is the relative change in the binary's binding energy 
caused by an encounter. This change, for an eccentric binary 
with a large mass ratio of one star to a companion and an 
encountering star, if occurr ed at the distance q(q^> a), can 
be ro ughly estimated as (e.g., lHeggiel il975HRoy & Haddowl 
l2003h 



mi I a 

mm \q 



3/4 



-(q/aY' 



(17) 



Assuming that all hardening happened with the maximum 
8 ~ 1113 /in^h ~ 0.04, it can therefore be estimated that, the 
hardening of a 1000 R Q BH-WD binary to, e.g., 35R Q should 
take about 100 encounters. As the collision time for a 1000 
Rq BH-WD binary is only about 10 7 yr, and is still only about 
10 8 yr for a 35 Rq BH-WD binary, it is plausible therefore that 
all 1000 R Q BH-WD binary can be hardened to 35fl Q within 
a few Gyr. We choose here 35Rq as the important separa- 
tion at which triple formation starts to have 100 % effi ciency 
in making this binary a MT BH-WD binary (see §2.31 . 

However, this is an idealistic picture. There are plenty of 
encounters that this binary would have undergone in order to 
reach 35 Rq, and certainly not all them will be simple fly-by 
hardening encounters. Some of the encounters would result 
in exchanges (where the lighter companion is often exchanged 
with a more massive incoming star, and the binary gets wider), 
mergers, or even binary ionizations, reducing therefore the 
chance for a binary to harden to the separation we are inter- 
ested in. 

We set up a numerical experiment, by taking a BH-WD 
binary with the specific initial binary separation a,. At the 
collision rate predicted by the current binary separation and 
the adopted n c = 10 5 pc" 3 , we bombard the binary with single 
stars drawn from the simplified core's population, randomiz- 
ing impact parameters and using Fewbody. Initial eccen- 
tricities of BH-WD binaries for this experiment is adopted 
to have a thermal distribution. A simplified core pop ula- 
tion is taken from Monte Carlo runs dlvanova et al. 2008|) and 
has only 4 different groups for the non-BH population: stars 
with masses 0.22,0.506,0.75, and 1.13M , with according 
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TABLE 1 
Hardening 





Hardened fraction [per cent] Average time [Gyr] 
Black hole fraction in the total core population [per cent] 
0.4 0.04 0.004 0.4 0.04 0.004 

Corresponding /bh.0.1 
10 1 0.1 10 1 0.1 






Hardening to 35 Rq 












Original binary survival only 






1OCU?0 


4.6 


12.9 14.8 14.9 0.86 


1.38 


1.48 


1.52 


25Oi?0 


0.53 


2.5 3.4 3.4 1.07 


1.9 


1.9 


1.9 


500S Q 


0. 


0.8 1.0 1.3 


1.9 


2.1 


2.3 






Exchanges are allowed 








lOOi?0 


9.4 


30.4 38.1 38.9 1.06 


1.50 


1.58 


1.60 


25Oi? 


1.3 


9.8 13.9 13.3 1.17 


1.85 


2.00 


2.07 


5OOi? 


0. 


3.2 4.8 5.9 


1.78 


2.15 


2.09 






Hardening to 80 Rq 












Original binary survival only 






100i? Q 


37.9 


56.0 58.2 60.6 0.22 


0.29 


0.30 


0.31 


15Oi? 


12.7 


27.5 30.0 30.9 0.35 


0.52 


0.57 


0.57 


2OOi? 


6.5 


16.9 19.1 20.7 0.41 


0.64 


0.68 


0.70 


25O# 


3.5 


11.3 14.4 14.4 0.46 


0.73 


0.77 


0. 77 


5OOi? 


1.1 


3.1 3.8 3.8 0.37 


0.90 


0.98 


0.98 



The table shows what fraction of BH-WD binaries with 15M0 BH and 
0.6MqWD and initial binary separations a\ can be hardened to 35 Rq and 
80 Rq, as a function of the relative BH content in the core: 0.4%, 0.04%, 
0.004% of the population, or no any other BHs present (/bh.0.1 are the 
corresponding values of the fraction of BHs that are available for interactions 
with normal stellar population, see §1). The average time corresponds to the 
time that it takes to harden those binaries that successfully got harder. The 
number of encounters used for the statistics: 10000 for 100^0, 5000 for 
100R Q and 1000 for 5OOR . 

33%, 38%, 25%, and 4% contributions to the non-BH popu- 
lation. Then we vary the BH component, starting from the 
highest possible, where the number of BHs is 0.4% of all other 
stars in the core (this is about 100% of the initially formed BH 
population, see the discussion in §1, or /bh.o.i = 10), 0.04% 
(this corresponds to the case of /bh.o.i = 1), 0.004%, and no 
contribution at all (this correspond to the case when only one 
BH, which is in the considered BH-WD binary, is left). Here, 
one BH corresponds to 0.001% of the core population. 

We consider the fraction of BH-WD binaries that would 
harden to 35/?© within 10 Gyr as well as the average time 
it will take them. We also separately consider cases where 
we insisted that the original BH-WD binary would made it to 
35^q, or simply any binary containing initial BH would made 
it (allowing multiple exchanges). The latter case could corre- 
spond to the case when almost all the stellar non-BH popu- 
lation in the core are WDs (e.g., in some runs performed in 
iFregeau eTa l. 2009b, WDs fraction in the core reached 70%). 
The results are shown in the Table Q] 

There are several trends in the results, and all of them well 
expected theoretically. Indeed, as expected, the fraction of 
survived and hardened BH-WD binaries /h ai d is well below 1 . 
The fraction /h al[ j decreases with increasing a\ and increases 
when the BH population in the core drops, as other BHs con- 
tribute significantly into exchange reactions or ionisation. The 
average time that it takes to harden a binary agrees with the 
analytic estimate above and is about 1-2 Gyr. Even though 
/hard increases when the BH population drops, it does not in- 
crease by a factor comparable to the relative drop in the BH 
population. Accordingly, since the resulting formation rate 
of hardened BH-WD binaries (their production from initially 
wider binaries) l\ a i-d oc /BH,o.i/hai-d, the highest Thard occurs in 



an unlikely case when most of initially formed BHs are re- 
tained (/bro.i = 10). The fraction /h arc j for /bh.o.i < 1 does 
not vary much with /bh.o.i ■ 

We also study what fractions can be hardened to a = 80/?©, 
for our optimistic scenario when all BH-WD binaries with a = 
SORq make it to the MT via TIMT. In this case, significant 
fractions of wide binaries can be sucesfully hardened. 

3. FORMATION OF A BH-WD BINARY 

In the previous Section we considered what can happen to 
a BH-WD binary in a dense stellar system. Here we ana- 
lyze which kind of a BH-WD binary could be formed and, 
accordingly, which formation mechanism dominates in the 
formation of potential BH-WD X-ray binaries. We neglect 
the possibility that an unperturbed primordial binary with a 
BH c ould be formed and s urvive in a dense stellar system 
(e.g., iDowning et al.l (120091) has shown that the fraction of 
primordial BH binaries is well below 1% compared to all 
BH binaries). It has been discussed that there are two pro- 
cesses through which a BH could acquire a non-degenerate 
stellar companion: exchange interactions with (primordial) 
binaries and tidal captures (for a thorough discussion of 
the latter mechanism with a BH and references therein, see 
iKalogeraet all 12004 . Here, we do not consider tidal cap- 
tures as they are inefficient for interactions between a BH 
and a WD; it was suggested that these encounters would re- 
sult r ather in a tidal disrup tion of a WD or its nuclear igni- 
tion (iRosswog et al.l 120091) . We consider exchange interac- 
tions, when an acquired companion is a WD, as well as two 
other possible mechanisms: physical collision with red giants 
(RGs) and three-body binary formation. The latter mecha- 
nism is usually neglected from consideration, as for all non- 
degenerate stars it is assumed that such encounters would 
rather lead to a merger than a formation of a hard binary. 

3.1. Exchange encounters 

In a binary that has been newly formed via exchange, the 
post-encounter binary separation will be a post ~ a pie niBH/mc 
where m c is the lower mass companion replaced by a BH. The 
mass m c would normally not exceed that of an average star in 
the core, and accordingly the post-encounter binary separa- 
tion will be at least 25 times larger than pre-encounter binary 
separation. In exchange encounters, post-encounter eccentric- 
ities are distributed thermally with the mean e ~ 2/3. 

From Fig. [T] we see that only those binaries with pre- 
encounter binary separation a pre < Q.3Rq could create a bi- 
nary that will evolve to MT directly (though up to a pre « 
3.1Rq if the final eccentricity is as high as 0.99), and with 
fl pre S 1.4/?q - when a formed binary will evolve to MT 
through a consequent triple formation. In both cases, a pre is 
small enough to expect a merger if at least one companion in 
the pre-encounter binary was a MS star. 

Therefore, to create a BH-WD binary with a post < 35Rq, 
only encounters with WD-WD binaries are relevant. The frac- 
tion of WD-WD binaries is just a few percent of the total bi- 
nary popula tion. From the simulations set in llvanova et al.l 
(120081 120061) we find the fraction of the short period WD-WD 
binaries at 10 Gyr in their "standard" globular cluster (corre- 
sponds to our typical dense globular cluster) to be fi, « 0.8% 
(as the fraction of the all objects in the core). Adopting 
fl pre w IRq, we find that the formation rate of potential BH- 
WD X-ray binaries through exchanges coupled with TIMT is 
Texch i ~ 2 x 10" 3 per Gyr per BH. 



7 



Relatively wide hard BH-WD binaries (a ~ 80- 1000fl o ) 
can be formed through encounters with both MS -WD and 
WD-WD binaries, as pre-encounter binaries can be up to 20 
Rq (here, MS-WD binaries could be only these that have 
a > 15Rq). Again, using the same set of simulations, we find 
that such binaries have a relative fraction fi, = 1 .4% and an 
average a pre ps 22Rq. The resulting formation rate is fairly 
similar to the formation via encounters with WD-WD bina- 
ries, r eX ch,2 ~ 0.06 /hard ~ 2 x 10~ 3 per Gyr per BH. 

In our "optimistic" case, when in all binaries with a post < 
80^q TIMT works, the fraction of binaries playing role for 
direct exchanges (with separations a < 3.27?© = 80,t?q/25) is 
/b ~ 2%. When averaged over all binaries in several numer- 
ical simulations, /h ar d ~ 25% (for /bh.o.i = 1)- We find then 
Texch.i = 6 x 10" 3 and r exc h.2 = 1 .5 x 10~ 2 per Gyr per BH. 

3.2. Physical collisions 

Physical collisions between compact stars (neutron stars, 
NSs) and low-mass red-giants can lead to a formation of NS- 
WD X-ray binaries (llvanova et al.l2005bl) . which, during their 
persistent phase, are also known as ultracompact X-ray bina- 
ries (UCXBs). From observations of UCXBs in the Milky 
Way's galactic globular clusters, we know that UCXBs are 
formed at a high rate. The fraction of UCXBs among all 
LMX Bs is much lar ger in globular clusters than in the field 
dDeutsch et al.H2000t) . 

When natal kicks from Hobbs et al. (2005) are adopted for 
NSs formed via core collapse and reduced kicks are adopted 
for NS formed via electron capture supernova, the theoreti- 
cally predicted UCXB formation rate through physical colli- 
sions is consistent with the observations of gl obular clusters 
both in the Milky Way amd in external galaxies ( llvanova et al.l 
I2005bl) . 

Like in the case of UCXBs formation, it is plausible that 
BHs, if retained and not detached in a BHs subcluster from 
the rest of the stellar population, will experience physical col- 
lisions with R Gs. The rate of such co llisions, per BH, can be 
estimated as in Iva nova et al.l (I2005bl) : 



Tbhrg ~ 2TTGf RG fpm BH n c R rg v c 



(18) 



'0.1/rg/p 



R RC lOkm/s 



15M 10 5 pc" 3 R c , 



per Gyr 



where Rrg is the average radius of the red giants, f v = r p //?RG 
describes how close has to be an encounter in order to re- 
sult in the formation of a binary compact enou gh to start the 
MT. For NSs, f p ss 1.3 dLombardi et alj|2006l) . However for 
BHs the maximum periastron that leads to RGs Roche lobe 
overflow at the closest approach is larger, / p w 5. 6 The frac- 
tion /rg of stars at the RG stage is typically ~ 0.4% of non- 
degenerate single stars at the age 10-12 Gyr, with the min- 
imum stellar mass in the popu lation be i ng 0.0 8M^ and with 
the initial mass function as in iKroupal J2002b . Mass segre- 
gation of stars in the core res ults in a higher /rg- Analyz- 
ing the set of simulations from llvanova et all (120081) . we find 
that this fraction is about 0.8% of the overall core popula- 
ti on and Rrg = 3-7 at 1 Gyr; in the simulations presented 
in iFregeau et al.l d2009al) . the post-processed /rg ~ 0.8% as 

6 We will show in the paper in prep, where BH-RG enc ounters are modeled 
using the SPH code described in Lombardi et al. 1 2006), that this is indeed 
the case, and the collisions with / p < 5 lead to the formation of bound binary 
systems. 



well. In our relatively dense stellar cluster, with « c = 10 5 pc 3 
and velocity dispersion of about 10 km/s, a BH of 15 M Q 
therefore could have ~3 x 10~ 3 / p chance of a physical colli- 
sions during 1 Gyr. Therefore physical collisions can provide 
a significant contribution to overall BH-WD X-ray binary for- 
mation only if all the collisions leading to the formation of 
bound binaries, with all the periastrons, up to 5Rrg, lead to 
X-ray binary formation. 

Let us estimate what fraction of the formed bound bina- 
ries can successfully become BH-WD X-ray binaries. A sim- 
ple estimate for the outcome of a BH-RG collision can be 
done using an energy argument, as in the standard common 
envel ope consideration, using the c*ceA formalism dWebbinkl 
119841) . Specifically, the final semimajor axis aj is given by 



(H = Rrc 



^BH ^RG,c( 
2 m RG WlRG.e 



(19) 



where mRG.core and /WRG.env are respectively the core and en- 
velope masses of the RG, and Rrg is the RG radius. The 
paramter A introduced to characterize the donor envelope cen- 
tral concentration can be found directly from stellar models 
and is about 1 for low-mass RGs. The parameter ckce i s i n ~ 
traduced as a measure of the energy transfer efficiency from 
the orbital energy into envelope expansion, and is bound to 
be below 1. The energy balance in the eq.(u"9l assumes that 
the amount of the energy that is transferred from the orbital 
motion to the envelope expansion does not exceed the energy 
that is required to eject the envelope to infinity. The kinetic 
energy that the ejected gas has at infinity is neglected, as well 
as that two stars had initially at the infinity, as it is very small 
compared to the binding energy of the envelope. 

It can be seen, if q:ce = L that the final binary separation for 
a case of a BH collision with a low-mass RG will result in this 
case in a formation of a binary with «f > 1.3/?rg (this mini- 
mum value is for a case of an early subgiant, when the core 
mass is minimum, e.g., it is ~ 0.12M Q for O.9M RG with 
~ 2Rq radius and z = 0.005). Comparison with the parameter 
space on Fig[T] shows that, if a post-collision eccentricity is 
low, only encounters with a subgiant /?RG'«RG.coi-e/'«RG,env ;$ 
O.66/?0 could lead directly to a BH-WD X-ray binary, this 
corresponds to a subgiant with mRG.core ~ 0.17M Q (consider- 
ing a O.9M RG; its radius then is only 2.75/? Q ). The time 
this RG spends as a subgiant before it reaches mRCcore ~ 
O.17M is only 135 Myr, or ~ 25% of its RG lifetime. As- 
suming that the final separation after a physical collision will 
be not exceeding the predicted by simple energy conserva- 
tion, the formation rates are then ~2x 10~ 3 / p per BH per 
Gyr. This rate is an upper limit, as it also is being assumed 
that all the physical collisions that lead to the formation of 
bound systems will lead as well to the RG envelope stripping. 
From SPH simulations, we find that collisions with / p > 1 will 
lead to a formation of binaries with e < 0.8, making / p w 1 to 
be a border line between these collisions that lead to the for- 
mation of BH-WD binaries able to start MT in isolation, and 
these that could not (see also Fig.l). 7 This non-TIMT channel 
therefore implies a formation rate of ~ 4 x 10~ 4 per BH per 
Gyr. 

For collisions with / p > 1, TIMT is the main mechanism 
that brings a formed binary to the MT. We can estimate the 
maximum optimistic fraction. For that, we find what fraction 

7 Details of the hydrodynamical simulations of collisions between BH and 
RG will be presented in the paper in prep. 
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of collision would lead to the formation of BH-WD binaries 
with a < 8O/?0, assuming that, even if a physical collision it- 
self could not form a binary that will start MT in isolation, 
TIMT will help (accordingly, / p = 5). Using the eq.(IT9b. we 
find that the RG should have /?RG'«RG,core/' w RG.env ^ 9.6R e . 
This occurs when a 0.9M Q RG has ntRG.core ~ 0.3 1M Q and 
^rg ~ 17.8^q and corresponds to almost 95% RG lifetime 
(accordingly, and roughly, all collisions with RG can be 
counted towards TIMT influenced collisions). We conclude 
therefore that most optimistic formation rate via physical col- 
lisions, when TIMT is taken into account, is 1.5 x 10~ 2 per 
BH per Gyr. 

3.3. Three-body binary formation 
Three-body binary formation occurs if 

• three objects will meet within some vicinity of each 
other a v ; 

• in the case of such a meeting at least two of the objects 
will form a bound system; 

• the formation of this bound system is not altered by 
physical collisions or strong tidal interactions between 
the objects due to finite size effects. 

In the literature, usually only the first condition is 
meant when gener alised three-body binary f ormation rate 
is desc ribed (e.g., iBinney & Tremainelll987t llvanova et alJ 
l2005ah . There are however limitations that make such 
treatments inapplicable to our case. The formation rate in 
IBinney & Tremainel (11987b is derived for equal masses only 
and is derived to estimate only the formation of binaries which 
have the hardness ratio rj of their binding energy to the ki- 
netic energy of an av erage object in the core to be 1. In 
llvanova et alJ d2005al) . the formation rate is derived to take 
into account unequal masses and different energies; however, 
the treatment is applicable only when the resulting binaries 
are near the hard-soft boundary or a bit harder. 

In their derivati on of the form a tion rat es of binaries with 
different energies, llvanova et al.1 (12005 al) . approximate that 
the orbital separation in the formed binary is the same as 
the size of the vicinity a v where the t hree objects meet. 
As sh own by numerical experiments in lAarseth & Heggiel 
(119761) . this assumption is approximately satisfied at least for 
equal masses, where a(l + c)/a v « 1. 

In our study, we are most interested in the formation rates 
where objects have a fairly large mass ratio and have high en- 
ergies (rj > 100). Compared to r\ = 1 binaries, the decrease 
in the hard binary formation rate due to physical collisions 
and tidal effects could be significant. We therefore limit our- 
selves to the consideration of only degenerate objects meeting 
each other. On the other hand, the limitation of the consid- 
ered cases to only very hard binari es eliminates the necessity 
to consider the second condition: Aarseth & Heggiel (119761) 
showed that the probability to form a binary is strongly in- 
creasing as vicinity is decreasing, P oc if 2 , and, even in the 
case of rj = 2, 77% of three-body encounters resulted in bi- 
nary formation. They have also shown that the average eccen- 
tricity of formed binaries is a bit higher than in the thermal 
distribution: < e >= 0.77 for rj = 2 cases, though this was the 
largest energy considered and the average values were slowly 
decreasing with hardness. We may expect that the average ec- 
centricities could go as low as that of the thermal distribution. 



The rate at which three objects meet in the same vicinity 
a v is the product of the rate Ti(ci v ,m\,ma) for two objects to 
meet in this neighborhood and the probability P3 that during 
this event a third object will pass by. As previously, the two- 
body encounter rate T2(a v ,mi,m2) is standardly derived as a 
combination of geometric cross-section and gravitational fo- 
cusing: 

2 V l 

T2(a v ,m u m 2 ) = Tm 2 v 00 a v (l + - r ) , (20) 

v 

where v? = 2G(m 1 + m2)/a Y . Two objects spend in this vicinity 
about t v w 2a v /v p . 

The probability that a third object will be within the same 
vicinity, P3, is then usually found only assuming that a third 
object will geometrically sweep a certain volu me during r v 
dBinnev & Tremainell987l ; [lvanova et al.l2005al) . A "geomet- 
ric" cross-section however is good only for cases of a v large 
enough so v p < Voq. In the case when we are interested in 
the encounters occurring within a small vicinity only, gravita- 
tional focusing must be taken into account. In some sense, this 
implies that we are looking for a probability that, at the same 
time, two objects will be gravitationally focused by a BH. The 
probability that a third object will pass within a v during r v is 
then 

#j~27rn 3 <4-^-. (21) 
The three-body formation rate then is 



r 3 (a v ) ~ 2TT 2 n 2 n 3 a 5 v ^- ~ 2 5/2 7r 2 G 3/2 « 2 n3fl"v^«4H- ( 22 ) 

For very hard binaries, this rate is significantly larger than 
the corresponding rate derived in llvanova et alJ d2005al) and 
can be written as 

T 3 (fly) =il.6x 10 _16 per Gyr x (23) 
( «c V ( M 3 ' 5 / 10km/s \ 2 / m BH \ 3/2 
Uo 5 pc-V \Rq) I J \15McJ 

Therefore, even in a relatively dense stellar cluster, a BH 
of 15 Mq will have only about 2 x 10~ 9 three-body encoun- 
ters within its 100 Rq vicinity during 1 Gyr. The formation 
rate necessary to explain BH-WD X-ray binaries formation is 
consistent only with encounters within ~ SQQQRq, or rj = 10. 
We also note that in this case v p ^> v^, the assumption used 
for the derivation, is no longer valid. Even in a denser clus- 
ter (n c ~ 10 6 pc) with a somewhat smaller velocity dispersion, 
the rate given by the above estimate is still about an order of 
magnitude less than the observed formation rate. 

We conclude that unless a significant fraction of small- 
vicinity three-body encounters lead to a formation of binaries 
much smaller than the vicinity where they met (10-50 times 
smaller), three-body binary formation can not explain the ob- 
served rates of BH-WD X-ray binary formation. 

4. CONCLUSIONS 

We have discussed which BH-WD binaries can become X- 
ray binaries, how these binaries can be formed and at what 
rate. We have considered the following formation channels: 
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TABLE 2 
Formation rates 



Channel 


Conservative 


Optimistic 


EX + TIMT 

EX + HARD + TIMT 

PC 

PC +TIMT 


2 X 1(T 3 
2 X 1(T 3 
4x 10 -4 


6 x 1(T 3 
1.5 x 1(T 2 

1.5 x icr 2 



The table shows the formation rates of BH-WD MT binaries, per BH per 
Gyr. Here "EX + TIMT" is formation via an exchange encounter with the 
consequent TIMT (triple-induced mass transfer), "EX + HARD + TIMT" is 
an exchange with a wide binary formation with the consequent hardening and 
TIMT. "PC" is a physical collision directly leading to MT, and "PC + TIMT" 
is the formation through a physical collisions of a bound but wide BH-WD 
binary that is brought to MT by TIMT. "Conservative" is the case when we 
take into account only BH-WDs with a < 35i?©, which are guaranteed to 
create Kozai triple that results in TIMT during 1 Gyr. "Optimistic" is the 
case when all BH-WD binaries with a < 80Rq have a chance to be in such 
Kozai triple that leads to TIMT, a process that takes several Gyr. 

• exchange encounters 

• physical collisions with giants 

• three body binary formation 

All of the channels require that at least a fraction of BHs in- 
teracts strongly with other stars in the cluster. 

We find that the most important mechanism to make a 
BH-WD X-ray binary from an initially dynamically formed 
BH-WD binary is the triple induced mass transfer (TIMT), 
the second most important is hardening with the consequent 
TIMT (see Table 2). TIMT is found to be very effective for 
BH-WD binaries and is not expected to be effective at all for 
BH-MS binaries, as most of potential triple formation encoun- 
ters will result in BH-MS collision. We do not find physical 
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With regard to TIMT, we distinguish between "conserva- 
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we take into account only BH-WD binaries with a < 35/?q. 
These binaries are essentially guaranteed to undergo TIMT 
within 1 Gyr. In the optimitic case, we take into account all 
BH-WD binaries with a < 80R Q . For these binaries, TIMT 
has a time-scale of several Gyr. 

With our conservative estimate of the formation rates, we 
find that we can explain the formation rate inferred from ob- 
servations, ~4x 10~ 3 //BH.o.ii but only if /bh.o.i ~ 1- This, 
by our definition of /bh.o.i, corresponds the case when 10% 
of all formed BHs remain in the cluster and interact with the 
core's stellar populations. 

With the optimistic scenario, we find that all the channels 
have comparable formation rates and each of the channels can 
explain the observed formation rate of BH-WD X-ray bina- 
ries. The combined rate of all the channels is ~ 4 x 10~ 2 per 
BH per Gyr. Comparing this with the formation rates inferred 
from the observations, we find that we can explain the ob- 
servations only with /bh ^0.1. In other words, even in the 
most optimistic case, we require that at least 1% of all ini- 
tially formed BHs (/bh.o.i =0.1) should be both not evapo- 
rated from the cluser and not dynamically detached from the 
core's stellar population to BH subcluster. Future simulations 
could address this. 
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